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ABSTRACT 


Aims. In this study, we conduct a pilot program aimed at the red supergiant population of the Magellanic Clouds. We intend to extend 
the current known sample to the unexplored low end of the brightness distribution of these stars, building a more representative dataset 
with which to extrapolate their behaviour to other Galactic and extra-galactic environments. 

Methods. We select candidates using only near infrared photometry, and with medium resolution multi-object spectroscopy, we 
perform spectral classification and derive their line-of-sight velocities, confirming the nature of the candidates and their membership 
to the clouds. 

Results. Around two hundred new RSGs have been detected, hinting at a yet to be observed large population. Using near and 
mid infrared photometry we study the brightness distribution of these stars, the onset of mass-loss and the effect of dust in their 
atmospheres. Based on this sample, new a priori classification criteria are investigated, combining mid and near infrared photometry 
to improve the observational efficiency of similar programs as this. 


1. Introduction 

When stars with masses between ~ 8 and ~ 25 Mq deplete the 
hydrogen in their cores, they quickly evolve away from the main 
sequence to the cool side of the Hertzsprung-Russell diagram, 
becoming red supergiants (RSGs). Evolutionary models indicate 
that this change happens at almost constant bolometric luminos¬ 
ity. Therefore the decrease in temperature has to be compensated 
by a rapid and very significant expansion of the atmosphere, 
which reaches a radius in the range of 400 to 1500 Rq. 

The age range for stars in the RSG phase goes from 8 Myr 
to 20 Myr ( |Ekstrom et al.|2013] l. All RSGs are thus young stars, 
associated with regions of recent star formation. However, the 
nature of the initial mass function constrains the number of stars 
able to evolve into RSGs: [Clark et ah] ( |2009^ estimated that a 
cluster must have a minimum initial mass approaching 10"^ Mq 
to guarantee the presence of 2 -4 RSGs at a given time. 

The RSG phase is critical to our understanding of high-mass 
star evolution. This phase lasts < 10% of the lifetime of a star, 
but the physical conditions, specifically mass loss, as a RSG will 
determine its final evolution. In consequence, evolutionary mod¬ 
els for high-mass stars find a strict test-bed in the RSG phase. 
However, the physics of RSGs defies the limits of ID compu¬ 
tations, and much more complex models are needed to explain 
observations ( jEkstrom et al.|20f3] l. 

The interest of RSGs goes beyond their role as evolution¬ 
ary mod el constraints. Their high lumi nosity, from ^ 10"^^ to 


10^ ^ Lq ([Humphreys & Davidson||1979ll, combined with the fact 


that their emission peaks in the near infrared (NIR) allows their 
observation at these wavelengths out to very large distances, 
even if they are affected by high extinction. Thanks to this, in the 
past few years several massive and highly reddened clusters have 
been discovered in the inner Galaxy (Clark et al.|2009bl Davies 


jet al.||2007| |2008| |2012[ [Negueruela et al.||2010| ). As their only 

observable components are the RSGs, these clusters are known 
as red supergiant clusters (RSGCs). 


touching the base of the Scutum-Cmx arm ( [Davies et al.|200'7] l. 
This region contains at least three large RSGCs, but as these 
works show, there are many other RSGs around these clusters, all 
with similar radial velocities, suggesting that they all belong to 
the same dynamical structure. Eor this search, many candidates 
were selected based on photometric criteria, but other late-type 
luminous stars, such as red giants or asymptotic giant branch 
(AGB) stars, have similar photometric characteristics. Therefore 
the only way to disentangle these populations is by looking at 
their spectra. 

The aim of the present study is to extend our search for new 
RSGs to other galaxies. Contrasting extragalactic samples with 
those of the Milky Way will provide information about the pro¬ 
cesses of stellar formation and evolution. Also, [Davies et al.j 
( 2010[ l showed that RSGs can be used as abundance probes, 
opening up a new method to study the metallicity in other galax¬ 
ies. Einally, the physical characteristics of RSGs vary for differ¬ 
ent metallicities ([Humphreys|1979a[ Elias et al.|1985[ Levesque 


[et al. [20061 [Levesque & Massey|2012[ l, but the number of known 

RSGs beyond the Galaxy and the Magellanic Clouds (MCs) is 
very low ( Levesque|2013[ l, and there is not a complete sample of 
RSGs in any Galaxy, not even the relatively nearby MCs. 

In this extragalactic search, our first step are the MCs, where 
large RSG population has already been studied ([Humphrey 


1979a|b[ Prevot et al. 1983[ Elias et al. 1985[ Levesque et al. 


2006 1 . This has the big advantage that the distances to both 


Recently, Negueruela et al. (2011 2012) searched for new 
RSGs in the region where the end of the long galactic bar is 


clouds are well established, removing part of the degeneracies 
that plague the studies of RSGs in the Milky Way. However, un¬ 
til now, only the brightest RSGs of the MCs have been studied, 
leaving almost unexplored the dimmer end of the RSG popula¬ 
tion and the frontier between RSGs and AGBs. 

In the last fifty years many works have studied photomet¬ 
rically the red population of the MCs, the high-mass popula¬ 
tion, and the RSGs themselves. The most recent and exhaustive 
among these works are Massey| ( [2002 1 and Bonanos et al. ( 2009| 
[2010[ l. Massey ( 2002) surveyed the MCs in the visible, looking 
for high-mass stars. He found ~280 RSGs in the Large Mag- 
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ellanic Cloud (LMC) and ~160 in the Small Magellanic Cloud 
(SMC), many of them previously unknown. However, as RSGs 
have their emission peak in the near infrared and their high mass- 
loss tends to redden them, by using visible data only the the less 
reddened and/or optically-brightest RSGs were observed. 

Bonanos et al. explored the mid infrared (MIR) properties 
of high-mass stars in the LMC ( Bonanos et al.|2009] l and SMC 
( |Bonanos et al.|20T0| l. Using spectral classifications taken from 
the literature, they derived new MIR criteria to identify RSGs. 
Britavskiy et al.|(20I4[) used these criteria to select a few dozen 


candidates in other galaxies, finding six new RSGs. 

There have also been some spectroscopic surveys aimed at 
these populations. The first works (|Humphreys||1979¥ 


_ _^ pr 

et al.|1985]l studied a small number of the brightest RSGs, con 


firming their nature. More recent works have taken advantage 
of the availability of large scale photometric surveys to select 
large numbers of candidates that can be observed efficiently us¬ 
ing multi-object spectroscopy. Massey & Olsen] ( f2003[ ) obtained 
spectra for a statistically significant sample taken from their pre¬ 
vious photometric survey, confirming the RSG nature of most of 
their candidates. [Neugent et al.| ( |2012) l did a selection of RSGs 
and yellow supergiant (YSG) candidates using 2MASS. Of their 
1949 RSG candidates, they observed 522, labelling 505 of them 
as "probable supergiants", as even when their apparent magni¬ 
tudes and radial velocities were compatible with membership to 
the LMC, the authors did not perform any spectral classification. 

In this work, we present a larger sample of candidate RSGs 
in the MCs, all of which receive a spectral classification based 
on intermediate-resolution spectra. We analyse the success rate 
of different photometric selection criteria and discuss the pos¬ 
sibility of generating clean samples of RSGs. In addition, we 
take advantage of photometric catalogues to study the spectral 
energy distributions of all candidates and search for further se¬ 
lection criteria. Finally, we discuss statistical properties of the 
RSG population in the MCs. 


2. Target selection and observations 

2.1. Overall strategy 


Observations were carried out with AAOmega at the Anglo- 
Australian Telescope as backup/filler for a longer programme 
aimed at the inner disc of the Milky Way during 9 nights between 
2010 and 2013. The observational strategy behind this comple¬ 
mentary program was built upon three different samples: a set 
of photometrically selected candidates, a group of previously 
known RSGs (from |Elias et ar]|1985[ [Massey & 01sen| 


2003 


and a third group of known YSGs from Neugent et al. (2010 1 


used as low priority targets in areas of low target density (only 
in the SMC, as in the LMC we constrained ourselves to an area 
of high target density). 

The samples of already known supergiants (SGs) were added 
partly as a control sample to compare with our new candidates 
and partly to extend classification criteria and schemes from the 
blue and optical into the Can triplet range at low metallicity. 
This paper deals mainly with the new candidates, while the full 
potential of the other samples will be exploited in future publi¬ 
cations. 


2.2. Selection criteria 

The recent boom in the detection of RSG rich clusters is mainly 
due to the availability of all-sky, near infrared data, because in 
this wavelength regime the peculiar extended atmospheres of 
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Fig. 1. Proper motion (taken from USNO-Bl) relative distribution for 
putative RSGs of the SMC (black line) and for all the stars in the field 
(red line). 


these stars stand out. In fact, just by using 2MASS photometry 
( Skrutskie et al.|[2006 1, it is possible to define a pseudo-colour 
Q = (J - H) - l.S X (H - Ks) able to separate between blue and 
red stars. This parameter has been proved to be a excellent tool to 
pick out RSGs, as they often show values of Q similar to those of 
yellow stars (Q ~ 0.2 - 0.3). This property, combined with their 
unusual brightness in the Ks band, allows the definition of purely 
photometric filters to select this population while minimizing in¬ 
terlopers. These criteria, combined with spectroscopic follow-up 
to confirm the nature of the candidates, have been used success¬ 
fully and extensively in our own galaxy ([Negueruela et al.|2()TT 


|2012[jGonzalez-Fernandez & Negueruela|2012| l and in this work 

we extend their use to the Magellanic clouds. 

The fact that these dwarf galaxies are not part of the Milky 
Way has the advantage that foreground objects are more eas¬ 
ily filtered out, particularly if they have large measured proper 
motions. Background objects, in contrast with the disc of our 
galaxy, are scarce and fall outside the parameter space occupied 
by RSGs. This comes at the price of a larger distance modulus, 
but as RSGs are intrinsically very bright, this is not an important 
issue. Also, the reddening towards the clouds is relatively small, 
with typical values around E{B-V) ~ 0.1 ( [Soszynski et al.|2002[ 
Keller & Wood|2006)l and so the pseudo-colour Q, that relies on 


the assumption of a given extinction ratio between bands, will 
not be affected by variable or non-standard extinction laws, at 
least outside the most reddened sites of recent stellar formation. 

With all these considerations in mind, we defined a set of 
selection criteria for RSG candidates as follows: 


- Candidates should have 0.1 < 2 < 0.4. 

- They have proper motions compatible with the MCs (see Fig. 

0 - 

- The brightness divide between RSGs and AGBs is not well 
established, but RSGs show normally absolute magnitudes 
brighter than -8 in the Ks band. Taking into account the 
distance modulus to the clouds, these stars should appear 
brighter than Ks - 11. 

- Lastly, to optimize the observing time, we impose a cut at 
OT/ = 13 so that spectra with high enough signal-to-noise can 
be obtained in less than 30 minutes. 
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Fig. 2. Spatial distribution of targets in the SMC, over a DSS-Red image 
roughly 3° x 3° in size. 



Fig. 3. Spatial distribution of targets in the LMC, over a DSS-Red image 
roughly 4° x 2° in size. As can be seen, all the targets (from a single 
AAOmega pointing) are distributed over a region that covers less than 
50% of the galaxy. 


As a subset of targets were observed using several configurations 
and some spectra did not reach usable S/N, our sample amounts 
to a total of 617 individual objects for the SMC, and 314 for the 
LMC. 


2.3. Observations 


2.4. Data reduction 


While traditional spectral classification criteria for stars are nor¬ 
mally defined over the blue end of the optical range, many 
works have extended them for RSGs to the wavelength range 
around the infrared Can triplet, as it contains several atomic 
and molecular lines of physical interest. With the fibre-fed dual¬ 
beam AAOmega spectrograph it is possible to cover both regions 
of the spectrum for several hundred objects in a single exposure, 
making it an ideal instrument for this kind of studies. As it sits on 
the 3.9 m Anglo-Australian Telescope (AAT) at the Australian 
Astronomical Observatory, it has good access to the low latitude 
fields of the Clouds while offering a collector area large enough 
to observe their RSGs in a very efficient manner. 

The instrument is operated using the Two Degree Field 
("2dF") multi-object system as front end, allowing the simul¬ 
taneous acquisition of spectra through 392 fibres. These fibres 
have a projected diameter of 2'.'1 on the sky and are fed into the 
two arms via a dichroic beam splitter with crossover at 5 700A. 
Each arm of the AAOmega system is equipped with a 2kx4k 
E2V CCD detector (the red arm CCD is a low-fringing type) 
and an AA02 CCD controller. While the red arm was always 
equipped with the 1700D grating, the blue arm changed between 
the 580V and 1500V gratings. A summary of the configurations 
is offered in Table [T] However, since the projection of the spec¬ 
trum from each fibre on the CCD depends on its position on the 
plate, it is not possible to give a precise common range for each 
configuration. This effect displaces the spectral range limits by 
~20 A in the red range, ~40A for the 580V grating in the blue 
range, and ~20 A for the 1500V grating in the blue range. 

The nominal resolving powers (A/6A) at blaze wavelength 
for the 580V and 1500V gratings are 1 300 and 3 700, while the 
1700D grating reaches R ~ 10000 around the Ca triplet, allow¬ 
ing the measurement of line-of-sigh velocities with enough pre¬ 
cision for our purposes. 

The main body of the SMC was covered with two pointings 
(Fig .[^ that were observed using 8 different configurations, for a 
total of 1448 spectra. Only one pointing was devoted to the LMC 
(Fig .[^, visited with two configurations for a total of 464 spectra. 


Data reduction was performed using the standard automatic re¬ 
duction pipeline 2d£dr as provided by the AAT at the time. 
Wavelength calibration was achieved with the observation of arc 
lamp spectra immediately before each target field. The lamps 
provide a complex spectrum of He-tCuAr-i-EeAr-i-ThAr-i-CuNe. 
The arc line lists were revised, and only those lines actually de¬ 
tected were given as input for 2d£dr. This resulted in very good 
wavelength solutions, with rms always <0.1 pixels. 

Sky subtraction was carried out by means of a mean sky 
spectrum, obtained by averaging the spectra of 30 fibres located 
at known blank locations. The sky lines in each spectrum were 
evaluated and used to scale the mean sky spectrum before sub¬ 
traction. 


2.5. Measuring v\os 

We measured velocities along the line of sight by calculating the 
correlation function of our observed spectra with a suitable tem¬ 
plate of known non-cosmological redshift. Eor late-type stars as 
the ones here studied, a high resolution spectrum of Arcturus 
is normally used, but while for the part of the spectrum with 
/I > 1 fim this is an adequate standard, as the overall shape of 
the spectrum does not change dramatically, this is not the case 
for the wavelength range around the Calcium triplet, as in this 
region the spectrum of successive populations will be dominated 
by the Paschen series, then the Calcium atomic lines and lastly 
TiO molecular bands. This results in a rather dramatic change in 
typology, making it very difficult to find a one-size-fits-all stan¬ 
dard to use as comparison. 

We have chosen instead to use a whole family of MARCS 
synthetic spectra taken from the POLLUX database ( |Palacios| 
|et al.|201^ . Eor each observed spectrum, the most similar model 
is selected doing a first pass over the whole set of synthetic spec¬ 
tra using very rough increments in velocity (Ayios = 10 kms^') 
and once the best match is selected, a refined value of yios is mea¬ 
sured calculating the correlation between observation and model 
using 0.3 kms“^ increments. These measured velocities where 


Article number, page 3 of 


39 







A&A proofs: manuscript no. gonzalezfemandez 


Table 1. Summary of the observations 


Year 

Nights 

Grating 

Blue arm 

-^cen (^) 

Range (A) 

Grating 

Red arm 

^cen (-'^) 

Range (A) 

2010 

3 

580V 

4500 

-2100 

1700D 

8600 

-500 

2011 

2 

1500V 

4400 

-800 

1700D 

8700 

-500 

2012 

1 

1500V 

5200 

-800 

1700D 

8700 

-500 

2012 

2 

580V 

4800 

-2100 

1700D 

8700 

-500 

2013 

1 

580V 

4800 

-2100 

1700D 

8700 

-500 



CT (km/s) 


Fig. 4. Histogram of the standard deviation for the measured Uios of re¬ 
peated observations. 


later transformed into the heliocentric system of reference using 
the rvcorrect package from Iraf. 

Using stars with repeated observations, we can obtain an es¬ 
timate of the total uncertainty in wios, including measuring errors, 
wavel eng th calibration, astrophysical noise, etc. As can be seen 
in Fig.|^ the typical velocity dispersion is around 1.0 km s', and 
we can assume a conservative 99% conhdence interval for our 
measurements of uios at 4 kms“*. Another source of dispersion 
that needs to be taken into account are possible systematic effects 
between different observing runs. Using all the available stars 
(main program and SMC/LMC backup) we checked for these, 
as we have repeated measurements for several objects. System¬ 
atic differences in Uhei were all below 1 km s * and have not been 
corrected, as some of the helds (particularly in the LMC) have 
very low redundancy and hence is difhcult to measure and cor¬ 
rect properly for this effect. 

In this article we will only use velocities in a relative sense, 
to discriminate between populations from the MCs and from dif¬ 
ferent Galactic components. Being so, we only worry about the 
internal consistency of the calibration, without the need of an an¬ 
chor point to check for systematics. In any case, as can be seen 
in Fig. 1^ the derived systemic velocity for the SMC is in very 
good agreement with the values from the literature, hinting at a 
very low systematic residual, if any. This is not the case for the 
LMC. Since we are not surveying the totality of the galaxy, we 
are heavily biased by its internal dynamic structure, and cannot 
readily compare with an "average" galactic velocity. 
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3. Results 


3.1. Spectral classification 

To classify the observed stars we used spectra obtained with the 
blue arm, as they cover the wavelength range where classical 
classification criteria are defined ([Morgan et al.||1943[ |Fitch &| 


Morgan|1951t Keenan & McNeil|1976[|Keenan & Wilson 1977 


Keenan & Pitts|1980| [Morgan et al.|1981||Tumshek et al. 1985 


Keenan|1987|l. These criteria were complemented with our 

secondary indicators, whose variation with spectral type (SpT) 
and luminosity class (LC) we derived through visual compari¬ 
son between the spectra of known standards. These were taken 
from the Indo-US spectral library (Valdes et al.||20041l and the 


MILES star catalogue ( |Sanchez-Bl^quez et al.||2006| l, and de¬ 
graded to our spectral resolution (roughly a FWITM of 2.1 A) 
when needed. 

[Humphrey ( |1979a| l reported that the metallicity differences 
between the Milky Way (MW), the LMC and the SMC do not 
change the behaviour of the atomic line ratios and other spectral 
features used in the spectral classihcation of RSGs. Therefore, 
it is possible to use MW standards as comparison, and the same 
criteria developed for RSGs in one galaxy are applicable to the 
others, as long as they are based on line (or band) ratios and 
not line (or band) strengths. With our extended sample, we can 
conhrm that there are no apparent differences in the spectra of 
RSGs from both clouds, even considering that our sample cov¬ 
ers a rather broad spread in spectral types. In consequence, we 
adopted the same criteria for both MCs, using Galactic standards 
as comparison. 

We have performed our own classihcation for all the stars 
observed, even those with early SpTs (most of them part of the 
YSG control sample). As this work centres around RSGs, we 
will only discuss the detailed classihcation of stars with spectral 
type later than GO. This also avoids the metallicity effects over 
the classihcation of earlier spectral types, that is more heavily 
affected by this parameter (cf. |Evans & Howarth|2003) l. We have 
also found some carbon and S-stars among our targets, but these 
are easily identihed due to their very characteristic spectra. As 
these are interlopers in our sample of new candidates, we did not 
perform any further analysis on them. 

Spectra observed with the 580V grating cover roughly from 
3730 A to 5850 A (the exact limits depend on hbre position), 
but the S/N blueward of ~ 4500 A is very low for many of our 
stars (this is not the case for the earlier spectral types). In con¬ 
sequence, most of our classihcation criteria lie between 4500 A 
and 5850 A. 

The main LC indicators we used are the ratios betwe en the 
the lines of the Mgi triplet (5167 A, 5172 and 5184 A) (Fitch 
|& Morgan] 195 1) 1. From GO to ~M3, Mg i 5167 is clearly deeper 
than the other lines for LC 1. These ratios change slowly with 
SpT, but this variation is not enough to complicate the identihca- 
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tion of mid- and high-luminosity SGs (lab, la). There are other 
spectral features that can be used to confirm the LC: the ratio of 
Fej+Y II blend at 4376 A to Fe i at 4383 A (Keenan & McNeil 


1976 I, the ratio between the blended Fe i lines around 5250 A 


and the Ca n-Fe n-Ti i blend at 5270 A ( [Fitch & Morgan 
and the ratio of Mn i 5433 A to Mn i 5447 A. 


1951 


As the Balmer lines decrease in strength with SpT while at 
the same time metallic lines become more intense, to identify G 
or earlier subtypes we used the ratio of the H/? and Hy transitions 
to other nearby metallic lines. For those stars with enough S/N in 
this region, we also compared Hy to the CH G-band (from 4290 
to 4314 A): F stars have deeper Hy, at GO Hy and the G-band 
have similar depths and from GO to mid G, the G-band becomes 
dominant. 


For later SpT, the shape of the metallic lines changes due 
to the appearance of TiO bands. The sequence of M types is 
defined by the depth of these bands ( Turnshek et al.|1985| l, but 
their effects are noticeable from KO onwards. Starting at Kl, a 
TiO band rises at 5167 A, changing the shape of the continuum 
between Mg lines. We used this change to obtain the spectral 
type for stars between Kl and M2. However, at ~M3, the band is 
so deep that the Mg lines are not useful any more. This is also the 
case for the other indicators mentioned before, whenever there is 
a band close to them. 


The effect of the TiO band at 5447 A over two lines, 
Mni 5447 A and Fei 5455 A, can be used to obtain the sub- 
type of K stars, and the the atomic lines and molecular bands 
between 5700 to 5800 A, to obtain the spectral type for mid and 
late M stars (Mn 1 5718 A, Mn 1 5718 A, V 1 5727 A, VO 5737 A, 
TiO 5759 A, Til 5762). 


For spectra observed with the 1500V grating, we had to re¬ 
sort to different criteria, but the methodology was the same. We 


identified th e LC using the following ratios from Keenan & Mc- 


Neil| ( |1976| ): Fe i+Sr ii 4216 A to Ca 1 4226 A, Fe n-Y ii 4374.5 A 
to Fe I 4383 A and Fe i 4404 A to Fe i-hV n-Ti ii 4409 A. In all 
cases, the ratio is ~1 for LC I and »1 for less luminous stars 
(LCIII-V). 


SpT can be evaluated by comparing Hd at 4102 A and Hy 
at 4341 A with nearby metallic lines. For early-G types the 
depth of Hy is similar to that of the G band, while for F or 


earlier types, Hy is clearly dominant. Even if the Fe i 4347 A, 
Fe n-Cr n-Ti ii 4344 A and Mg n-Cr n-Fe 1 4351 A lines vary with 
LC, they can be used to determine SpT, because we can constrain 
this parameter with other indicators. 

We also used the ratios between Fe 1 4251 A and Fe 1 4254 A, 
Fe I 4280 A and Fe n-Ti i 4282 A, and the behaviour of the lines 
Fei-nCoi 4579.5 A, Fei 4583 A and Fen-Cru-Cai 4586 A to 
confirm the spectral subtypes in the later G and K sequences. 

Despite the fact that the first TiO bands appear already for 
K type stars, it is only in redder spectral regions. In the spec¬ 
tral range considered here they are not noticeable until early M 
subtypes, and no TiO bands are clearly visible before ~M3. To 
classify the early M stars we used the spectral range from 4580 
to 4590 A and from 4710 to 4720 A. 


In order to facilitate calculations, we parametrized SpT and 
LC over a linear scale, assigning integers to each type and class. 
In those cases in which we doubt between two consecutive clas¬ 
sifications, we assigned the intermediate half-integer value. 

Even though the term RSG is normally applied to SGs with 
types K or later, we have also included G stars for two reasons; 
firstly, RSGs are intrinsically variable and some can change their 
type from early K to late G. Secondly, at lower metallicities. 


Table 2. Repeatability of our spectral and luminosity classification us¬ 
ing stars with multiple observations. 



A(LC) 

A{SpT) 

Number of stars 

2010 

0.4 

1.2 

99 

2011 

0.4 

0.8 

101 

2012 

0.4 

0.9 

129 

Weighted average 

0.4 

1.0 

329 


the typical SpT of a RSG becomes earlier. In consequence, as 
Levesque (2013|l noted, if we exclude G stars we are losing part 


of our target population (evolved high mass stars), specially for 
low metallicity galaxies such as the SMC. Therefore we have 
used all the stars with SpT GO or later for subsequent calcula¬ 
tions. 

There is some overlap between different SMC observations. 
As we performed the classification for each of the spectra of 
these redundant targets independently, we can use use them to 
test the internal coherence of our classification scheme. The final 
SpT and LC for these stars were obtained by averaging and using 
the S/N of each spectrum as weight, rounding the final figure to 
the closest entire or semi-entire number. 

The mean differences between the spectral classifications of 
these repeated targets are given in Table [^ As can be seen, the 
differences in both LC and SpT are of the order of the classi¬ 
fication step, as long as we take into account that we assigned 
semi-entire SpT only in those cases where the classification be¬ 
tween two consecutive subtypes was not clear. 

Attending to the obtained differences in our classification, 
we have assumed an uncertainty of +1 in SpT and +0.5 in LC 
for all our stars, even if there are no repeated observations for the 
LMC, as the observing conditions and the classification scheme 
were the same for all fields. 

Of all the stars with more than one observation, there are a 
few that present large discrepancies between epochs. Even if the 
numbers are compatible with normally distributed errors, we re¬ 
vised all these spectra to check the source of these differences. 
In many cases it is due to one of the spectra having low S/N. In 
these cases, as our final classification was done using the S/N as 
weight, the final result will be dictated by the high S/N classifi¬ 
cation. Other stars have good S/N in all their spectra, and differ¬ 
ences arise due to the lack of enough standards for some spectral 
subtypes and luminosity classes. This is the case of many G stars. 

A high fraction of RSGs are known to be long period vari¬ 
ables ( [Wood et al. 19831. These variations reflect not only on 
their brightness, but also on their spectra and radial velocities. 
Therefore, among many of our repeated observations, stars will 
appear with different SpT, LC and radial velocities from epoch to 
epoch. These classifications have also been averaged. As much 
of what is discussed in later sections deals with single epoch 
photometry, there is no use on retaining several classifications 
for the same object, and by averaging we ensure that the final 
values for SpT and LC will not be at any of the extremes, in¬ 
creasing the chances of better agreement with the asynchronous 
photometric measurements. 

However, for all those stars showing variability, we retain the 
different classifications, while for non-variable objects a single 
value is listed. We have considered to be variable all those stars 
which have a difference in SpT or LC, non attributable to other 
factors, larger than twice the uncertainty interval: 2 subtypes or 
1 luminosity subclass. 
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Fig. 5. Observed heliocentric velocities for all the sources in the LMC 
(red) and SMC (black). Over-plotted with dashed lines are Gaussian fits 
to the distributions, and their p is marked with a vertical dashed lines. 
For comparison, the blue vertical lines denote the systemic velocities 
for the clouds (taken from|Massey & 01sen|2003|l. 


Fig. 6. Plot of the apparent magnitude versus the Uhei for the observed 
sources. Squares mark candidates for the LMC, circles do so for the 
SMC and stars are left for Carbon stars. The colour coding denotes LC, 
see Sect.|3.1|for an explanation of the chosen parametrization. 


Finally, we want to stress that our spectral classifications are 
merely morphological. When we classify an object as a super¬ 
giant, we are simply stating that several significant spectral fea¬ 
tures in its spectrum look more like those of supergiant standards 
than those of giant standards. We are not making any assumption 
about the physics of the stellar interior. Even though there is a 
generally excellent correlation between spectral type and physi¬ 
cal characteristics, this does not always have to be the case. For 
example, recently |Moravveji et al.| ( |2013] l have presented evi¬ 
dence that a Her, an M5 Ib-II MK standard and anchor point of 
the classification system (because it is the high-luminosity stan¬ 
dard with a later spectral type), is an AGB star of only ~ 3 Mq. 


3.2. Membership to the clouds 


The velocity distribution of our potential SMC sources can be 
accurately modelled by a Gaussian distribution with parameters 
{p = 149.6 km s“', cr = 23.7 km s“'), while for the LMC these 
become (p - 271.4 km s"', cr = 15.3 km s“'), as can be seen 
in Fig. 1^ Based on this, an initial clean-up of the sample can 
be obtained using hard cuts in yhei, using +3cr as threshold. Yet 
the populations from the MW and the MCs cannot be separated 
based purely on dynamical criteria, as there are halo stars that 
show velocities compatibles with those of the clouds. This can 
be seen in Fig. where the sources are colour labelled accord¬ 
ing to their LC. Both for the LMC and the SMC there are stars 
of classes III to V within the dynamical envelope of the clouds 
but with apparent magnitudes incompatible with their distance 
modulus and spectral classification. The only way to weed out 
these interlopers is through detailed spectral tagging. 

On top of these MW populations, the transition from RSG to 
AGB is smooth, and so both very luminous AGBs and Carbon 
stars will appear in photometrically selected samples. Although 
these will indeed be part of the clouds, in order to ensure a pure 
sample of SGs, it is again mandatory to perform a good spectral 
characterization of the sources. Using both wios and the spectral 
classification, we can perform the last cleansing of the sample 
in order to produce a catalogue of SGs in the MCs. The results 
from the different stages of this process are shown in Table 
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The final sample contains a total of 160 SGs in the SMC and 
123 in the LMC. Of these, 70% are previously unknown SGs. 


4. Discussion 


4. 1. Selection efficiency and interlopers 


In Sect. 3.2 we show that only the combination of nios and spec¬ 
tral classification can separate sources from the MW and the 
MCs. Once this filtering has been done, we can proceed to anal¬ 
yse the efficiency of our selection criteria. 

As can be seen in Table of the 585 new candidates, 48% 
turned out to have LC Ib-II or brighter, and hence can be clas¬ 
sified as SGs. The ratio of success for our selection is of 53% 
for the LMC and 45% for the SMC. This small difference arises 
mostly due to the fact that we only covered the LMC with one 
configuration aimed at its main body, while for the SMC we also 
sampled the outskirts, where the relative density of bright inter¬ 
lopers is higher. 

Most of these interlopers turn out to be MW disc popula¬ 
tion, along with some high velocity halo stars. The majority of 
these could be removed by the application of a cut on proper 
motions based on a catalogue with more precise measurements 
than USNO.B-1. Proper motions from Gaia will allow samples 
almost clean of MW populations. Among MC interlopers, car¬ 
bon stars are particularly conspicuous. At low magnitudes, they 
are the main contaminants in the LMC. Some of these stars may 
be easily filtered out by looking at their colours, as for the MCs 
they reach (J - K^) > 2, while in our sample no SG is red¬ 
der than (J - Ks) < 1.6. But since typically, carbon stars will 
have (J - Ks) > 1.4 ( |Cole & Weinber^|2002[ ), there is some 
overlap between both populations, and much more so for fields 
under heavier extinction, where the SG population will be dis¬ 
placed into the red. This is the case too for mid-infrared colours, 
where the overlap is even more complete, and from having simi¬ 
lar colours, it follows that carbon stars and RSGs will also show 
similar values of Q. Being so, it is expected that these stars will 
appear in any survey aimed at RSGs with enough depth to reach 
the low luminosity end of the Ib population, as can be seen in 
Fig-IZl the fraction of recovered SGs drops below Mk = -9, 
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Table 3. Filtering of the original sample of candidates according to several criteria. In parenthesis are the carbon stars with Ui,ei not compatible 
with the clouds and the number of previously undetected SGs. 


Cloud 

Total 

uios filter 

LC filter 

Carbon 

Final sample 

LMC 

237 

203 

125 

48 (3) 

123 (70) 

SMC 

400 

179 

162 

10(0) 

160(128) 



Mk 


Fig. 7. Fraction of the final sample for both clouds made up by SGs (LC 
Ib-II or more luminous), giants (LC II or less luminous) and carbon 
stars. Distance moduli are taken to be 18.48 for the LMC and 18.99 for 
the SMC, from | Walker| j20 1 2] ! and |Graczyk et al.|j2014[ l 


reaching a point at Mk -8 where in fact interlopers dominate 

the sample. 


4.2. Completeness of the sample 


An item of utmost importance when considering a sample is that 
of completeness. In our case, this is delimited by the most re¬ 
strictive of the criteria outlined in Sect. |2.2[ the cuts in Q. As has 
been shown, introducing hard thresholds in Q leads to a low pro¬ 
portion of interlopers in the sample, but at the same time there is 
the chance it might leave out some SGs too. We can check this 
by comparing our a priori photometric selection with other sim¬ 
ilar programs. Of those available, the most complete and deep is 
|Boyer et al.| ( |2011) , with the added advantage that their classifi¬ 
cation is based on the NIR and MIR photometry of their sources. 


The results of our filtering applied to the objects from Boyer 
|et al.| ( 20TT| l can be seen in Fig.|^ At face value, this plot seems 
to indicate that we are missing a large fraction of the potential 


RSGs: of more than 3000 putative candidates flagged by Boyer 
|et al.| ( |2011| l in the SMC, we have only observed around 200. 
But there are a couple of points that we have to take into ac¬ 
count; firstly, stars labelled as RSGs in Boyer et al. ( 201 l| l extend 

down to otks = 12.5, a magnitude that translates to Mk^ -6.5 

at the distance of the SMC, hardly compatible with what is ex¬ 
pected for this population (see, for example, Elias et aT]|1985[); 
secondly, the spatial extension of the study conducted by Boyer| 
|et al.| ( |20TT| l is much wider than ours. If we take this into account 
and we impose that RSGs must have M^s < -8, the equivalent 
sample from|Boyer et al.|(|201 l|l is trimmed down to 479 candi¬ 


• • RSG 
■< . AGB 



-0.2 0.0 0.2 0.4 0.6 

Q 


Fig. 8. Q distribution of the objects from |Boyer et al.| ( j201 l| l, cover¬ 
ing only the SMC. The black dashed lines mark the boundaries of our 
filtering scheme. 


dates, of which around 200 have 0.1 < 2 < 0.4. A similar result 
is obtained for other photometric surveys as |Massey| ( |2002[ ): of 
the 288 candidates that have 2MASS photometry, 58% pass the 
Q filter; and also for Groenewegen et al. ( 2009| l, as 40% of their 
RSGs (stars with Mboi < -8.0) clear our cuts. At the other end 
of completeness, from the 21 RSGs identified by |Buchanan et al.] 
P006| l from a pool of objects with colours and fluxes at 8 pm 
coherent with evolved, massive stars, 17 are picked out by our 
selection criteria. 


As it can be seen in Fig. cutting out sources with Q > 0.4 
does a good job at removing AGBs from the sample (and late 
type giants, not plotted there), but at the price of also filtering out 
a large fraction of RSGs. While in the disc of our galaxy, where 
this kind of filters are mostly used, AGB and giant contamination 
is a serious concern, and so missing on a significant number of 
RSGs can be an acceptable price, in the MCs we have the added 
value of being able to separate a large fraction of these lower 
luminosity stars just by looking at their apparent magnitudes, so 
we need to develop finer selection mechanisms. 


Beyond its completeness, we can also check if our filtering is 
biased towards or against given spectral types. Using the control 
sample, composed of SGs taken from Elias et al. ( 1985|),[Massey 


|& 01sen| ( 2003| l and |Neugent et al.| ( |2bl0| l, as these objects were 
selected disregarding their Q and cover a wide range of spectral 
types. 


As can be seen in Fig.|^the behaviour of Q with spectral type 
is relatively complex, related to the variation of the intrinsic NIR 
colours of SGs. This is shown in Fig.[^ the change in {H - Ks) 
with SpT is similar for giants and SGs (essentially, a temperature 
sequence), although the latter tend to appear redder. This is not 
the case for (J - H), where the behaviour is markedly different; 
this is probably related to the fact that at shorter wavelengths. 
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Fig. 9. Q values for the control sample. Top panel: Fraction of these that 
pass the Q threshold, marked with dashed lines. 


the spectral energy distribution of SGs is dominated by the com¬ 
bined effects of dust (the interplay between scattering emission 
and auto-absorption; [Smith et al.||2001] l, molecular absorption 
bands and opacity (that drops steadily from J to H). These 
effects dictate the colour of the star, with a weaker dependency 
on its temperature. 

This implies that our homogeneous filter in Q will have dif¬ 
ferent a priori completeness depending on SpT. This is detailed 
in the top panel of Fig. where we plot the fraction of SGs from 
the control sample that clear our hltering: while for mid-type G 
SGs this criterion works very well (even if the relative abundance 
of these objects is low), its efhciency drops as the spectral type 
sequence progresses; for late G and early type K SGs it can reach 
a rather low ~ 30% completeness. For later types the fraction of 
stars inside our boundaries increases more or less linearly with 
SpT, keeping over ~ 75%, although at the very end of the M se¬ 
quence our sample is too poor to draw any conclusion. This ties 
back with the results of Fig. as we carried our test on the SMC, 
where the spectral type distribution works against the efficiency 
of our Q hltering, the comparison with |Boyer et al.| ( |201 l| l works 
as a sort of worst case scenario. Both in the MW and in the LMC 
the fraction of M type supergiants is much larger, and in partic¬ 
ular in the disc of our galaxy the RSG distribution peaks around 
M2 dElias et al.|1985| l. 

This varying selection efhciency for different spectral types 
is of paramount importance when devising surveys for galaxies 
other than our own. It has been shown that the average spectral 
type of the RSG population depends on metallicity ([Humphreys 


1979a[ [Elias et ah 1985[ Levesque et al. 2006[ Levesque & 


Massey|2012|l andso as Z decreases, more RSGs will have ear¬ 


lier types, moving slowly into the region where selection com¬ 
pleteness is worse. Although these effects are very difficult to 
evaluate a priori just based on photometric data, one clear solu¬ 
tion is to just open the accepted Q range. As we can see in Fig.[^ 
this would include in the sample an increasingly large number of 
interlopers; while for the MCs it would be possible to weed them 
out, this is not the case in other fields, galactic or extragalactic, 
and hence the need to develop new strategies arises again. 

The influence of dust in the variation of Q is further sup¬ 
ported by looking at the MIR colours of these stars, that in WISE 
( Wright et al.|201()| ) show {W\ - WA) excesses indicative of dust 
emission. In fact, (W1 - W4) turns out to be a good indicator 


Fig. 10. Evolution of the NIR colours for our control sample as a func¬ 
tion of spectral type. Red solid lines mark the intrinsic colours for giants 
of the same type, taken from[Straizys & Lazauskaiti[j2009^. 


of spectral types in our sample, and it is possible to refine the 
selection scheme on Q using this property. As can be seen in 
Fig. 11 there is a linear relation between Q and {W1 - WA), and 


so we can define an estimation of Q based on MIR colours by 
Gmir = 0.39 - 0.06 X {W\ - WA). Using this, and whenever 
WISE photometry is available, we can derive a Qmir to compare 
with Q. Eor our sample, the standard deviation of Q—Qmir is 0.1, 
and so it is possible to define a threshold to attain a given com¬ 
pleteness. This selection scheme has an important caveat: most 
late giants will fall within the same range as SGs, and while this 
is not a problem for other galaxies, where simple cuts in mag¬ 
nitude can weed out these stars, in the MW this strategy is not 
feasible. 

Beyond checking our completeness, we can use these cross¬ 
matches with other surveys to evaluate their selection efficiency. 
This is summed up in Table In this table, we detail those ob¬ 
jects in common between the listed surveys and our sample of 
new RSGs, and whether we confirm their SG nature in the case 
of photometric surveys or we have discrepant classifications in 
case of spectroscopic surveys. Two effects mediate the numbers 
in this table: firstly, we have removed the control sample from 
this calculation, as these objects are selected a priori knowing 
their stellar nature. Secondly, the spatial overlap between our 
survey and those in the table is not complete, and this fact limits 
the number of common targets. 

4.3. Photometric properties of the sampie 


In Fig. 12 we plot a CMD of the confirmed SGs in this work. 
As can be seen, using the criteria outlined in Sect. 2.2 we obtain 
a set of candidates that, while overlapping with previous works 
at bright magnitudes, allow us to extend the search for SGs to 
low brightnesses, in a region of the CMD relatively unexplored. 
Some of the candidates have already been spectroscopically con¬ 
firmed as SGs (Table 1^, but for homogeneity reasons and due to 
the variable nature of the spectra of RSGs, they were left as can¬ 
didates and observed anyway. 

Beyond comparing with previous studies, a homogeneous 
sample this size can be exploited to study the photometric be¬ 
haviour of RSGs. To do so we combine the photometry coming 
from 2MASS and WISE (see Table|^for a summary of the wave- 
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Fig. 11. Pseudo-colour 2 as a function of MIR excess for our control 
sample. The solid red line corresponds to the best linear fit, while the 
dashed lines denote the 2cr threshold, that contains 95% of the SGs. 

Table 4. This table summarises the overlap between our sample of new 
candidates and other works. In the upper panel, we show the objects 
in common with other purely photometric surveys, and the number of 
them that we have confirmed spectroscopically as SGs. In the lower 
panel, we perform the comparison for spectroscopic surveys, indicating 
the number of SGs and other populations in common. In parenthesis are 
indicated those cases in which our classification contradicts that of the 
original work. 


Photometric surveys 


Paper Cloud Candidates Confirmed 



Westerlund et al. ( 

19811 

EMC 

39 

29 



Prevot et al 

■(W 

83 1 

SMC 

36 

33 




Massey ( 

i002 


SMC 

8 

4 




Massey ( 

2002 

1 

EMC 

6 

5 

Groenewegen et al. 

2009 

1 SMC 

7 

5 

Groenewegen e 

t at. 

(2009 

1 EMC 

4 

1 


Boyer et at 


11^ 

SMC 

108 

99 


Spectroscopic surveys 




Paper 



Cloud 

SGs 

Other 


Elias et al. 

1985 

1 

SMC 

18 

0 


Elias et at. 

1985 

1 

EMC 

7 

0 



OSK199S^ 



EMC 

3 

0 



Massey (2 

003 

1 


SMC 

0(1) 

7(6) 


Buchanan et a 

r (12006 1 

EMC 

0(1) 

0 


JNeugent et at 

7(5^ 

EMC 

28 

8 


References. (l) |Oestreicher & Schmidt-Kaler| ( |1998| > 


length coverage of these surveys). It is worth noting that these are 
both single epoch surveys, and a fraction of the stars on our sam¬ 
ple (AGBs, RSGs, etc.) are expected to be variable. This implies 
that not only photometry and spectral classification will be asyn¬ 
chronous, but also different bands have been observed at differ¬ 
ent epochs. Yet this variability is known to decrease in amplitude 
with wavelength, and as it is discussed in Robitaille et al. ( |2008| l 
the most extreme variables are reduced to amplitudes of a few 
tenths of a magnitude in the MIR. Even so, some dispersion due 
to this effect is expected in colour-colour and colour-magnitude 
diagrams. 


Fig. 12. NIR CMD of all the RSGs in the SMC, where our more com¬ 
plete spatial coverage allows for more meaningful comparison with pre¬ 
vious works. Blue dots mark those already present in the literature, 
while red dots do so for those confirmed in this work. The void around 
otk ~ 10.5 is an instrumental effect. 


Table 5. Details of the bandpasses used in the text. 


Eabel 

Survey 

disc (Mm) 

AT ifim) 

J 

2MASS 

1.235 

0.162 

H 

2MASS 

1.662 

0.251 

Ks 

2MASS 

2.159 

0.262 

[B^l] 

WISE 

3.353 

0.662 

[B^2] 

WISE 

4.603 

1.042 

[B^3] 

WISE 

11.561 

5.507 

[W4] 

WISE 

22.088 

4.101 


As has been mentioned, the magnitude of these RSGs in the 
J band is affected by several atmospheric effects, such as molec¬ 
ular opacity and the appearance of dust, that are controlled by 
the effective temperature and surface gravity of the star. Inter¬ 
estingly, [VFl] is expected to be mostly photospheric and not 
subject to a strong absorption by the outer layers of the stellar 
envelope. The interplay of these factors results in the fact that 
the (J- [VTl]) is a very good indicator of spectral type, as can be 
seen in Fig. There is an almost linear relation between SpT 
and this colour. There are clear hints that this behaviour saturates 
around M3, but our sample lacks stars of later type, so we cannot 
confirm this fact. 


4.4. Dust, and mass loss 


One of the most relevant physical phenomena affecting the at¬ 
mospheres of late-type stars is that of mass loss. [Josselin et al.| 
( |2000) l show over IRAS data that the (Ks - [12]) colour is a good 
measure of as mass loss, as this process will reflect on the MIR 
excess. The Vk3 band from WISE is the one that mimics most 
closely the IRAS 12/rm band, and so we explore the presence of 
mass loss using the (A’s-[VT3]) colour. As can be seen in Fig. 14 
all of the MW population present in our sample falls in a stripe 
of 0 < (Ks - [VT3]) < 0.3. Almost all SGs, AGBs and carbon 
stars present redder values of this colour. 

Although the bulk of carbon stars tend to be redder than SGs 
(the majority of which satisfies {Ks - [W3]) < 1.0), these stars 
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Fig. 13. Variation of (J - [Wl]) with the spectral type of the RSGs in 
our sample. 


never reach values beyond {Ks - [VF3]) = 2.0, while several SGs 
and AGBs do so. This probably hints to an enhanced mass loss 
and the presence of strong winds. Although the mass loss in SGs 
is supposed to be much stronger than in AGBs, the distribution 
of the latter in this plot seems to follow closely that of late-type 
SGs. As due to our selection criteria we are only picking the 
bright end of the AGB population, this similarity points again - 
even if with low statistical significance - to the fact that there is 
not a sharp-cut transition from AGBs to RSGs. Similar effects 
have been observed by other authors, e.g. |Yang & Jiang|P012| l. 

The number of SGs showing evidences of mass loss becomes 
important around K5, but even for the most evolved types there 
are still stars that show small MIR excesses (and so potentially 
low mass loss). At all types, bright luminosity classes tend to 
feature higher values for (Ks - [VT3]). 

The onset of dust in the outer layers of their atmospheres 
is one of the most relevant factors that dictate the photometric 
properties of cool, late type stars. Large granular compounds 
(complex carbon molecules, silicate grains, ice particles, etc.) 
form and coalesce in the outer layers of their extended atmo¬ 
spheres, and mediate their appearance in the NIR and MIR. 
The thermal emission from these particles becomes dominant 
at long wavelengths ( [Smith et al.|200l| ). In Fig. we compare 
([1T3] - [VT4]), a colour that should be a good indicator of the 
dust temperature in the outer layers of the stellar shroud, with 
(Ks - [1T3]), that measures mass loss. As noted before only late- 
type SGs show mass loss, with a marked tendency to higher val¬ 
ues than AGBs or carbon stars. 

As we can assign membership to the clouds to our stars, it 
is possible to estimate their intrinsic luminosities by assuming a 
given distance modulus to each cloud. We use p - 18.48 for the 
LMC ( |Walker|2012| ) and p = 18.99 for the SMC ( jGraczyk et ah] 
|2014| l. Although not all the stars will be at the same distance, 
the spread of both clouds in distance modulus is low, and the 
dispersion introduced in the intrinsic magnitudes negligible for 
the kind of qualitative analysis we are going to conduct here. 

In Fig. f^we plot two CMDs for our sample. 


As expected (Fig. 16 right panel), there is a brightness limit 

for SGs, and no star in our sample goes beyond Mj-10.5, a 

limit that translates into ~ -11.5. Although the brightest 
stars in these NIR bands tend to be of late type, there is not a 
strong correlation between brightness and spectral type. As ex- 
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Fig. 15. Colour-colour diagram for our sample of MC stars, representing 
Ik's - [Vf3]) (related to mass loss) against ([Vf3] - [Vf4]), related to 
the temperature of the outer dust layers. In this plot, symbol size is a 
function of LC (i.e larger symbols imply brighter classes). 


pected, AGBs and carbon stars occupy the lower end of the mag¬ 
nitude ladder, although some of the latter can reach Mj-9. 

This is not the case for the brightness in the MIR. SGs with 
types earlier than K all cluster around M[W 4 ] ~ -10. The K 
sequence of types is more or less evenly distributed between 
-10 > M[w 4 ] > -13. Finally, late type, M RSGs reach up to 

M[w 4 ]-15, with what appears to be a linear relation between 

absolute magnitude and colour, likely to arise from the increase 
of dust in their outer layer, as this would make the star dimmer 
in J (hence redder) while brighter in [W4], as this band is dom¬ 
inated by the thermal emission of dust. There are some stars of 
various types that break this upper limit in M[W 4 ], but the pho¬ 
tometry of all of them seems to be affected by the nebular emis¬ 
sion of large, nearby H ii regions, such as 30 Doradus. 

It is worth noting that both CMDs sample very different 
physical elements of the stars, as for objects with extended en¬ 
velopes and strong losses, in the NIR infrared we are looking at 
the central object while in the red part of the MIR these extended 
atmospheres are the dominant component. 

4.5. Objects of particular interest 

- HV 838; This SMC star is a known large-amplitude pho¬ 
tometric variable with a period of ~660 days ( [Wood et al. 
|1983| ). It has been considered both an RSG and an AGB 
star in different works. We have observed it in all three 
SMC campaigns but only the blue spectrum from 2011 has 
a SNR sufficient to allow a proper classification. The spec¬ 
trum shows strong H line emission. In the infrared region, it 
displays strong inverse P-Cygni profiles in the Ca triplet and 
other nearby lines. This behaviour is completely atypical in 
RSGs. 

In view of this, we decided to study the infrared spectra from 
2010 and 2012. There are no emission lines in these, but the 
spectral types are very late, explaining why the blue spectra 
have such a low SNR. Based only in the infrared spectra, we 
find types ~M7 II in 2010 and ~M8 II in 2012, but in 2011 
the infrarred spectrum reveals a late K Ib star. 

While this extreme change in spectral type is very unfre¬ 
quent among RSGs, jSoszyhski et al.j (|201 l|l found a varia- 
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Fig. 14. Left:( J - H) vs. (A's - [1F3]) diagram for the whole sample. Right: (A's - [VF3]) as a function of spectral type for all the SGs in the sample. 



Fig. 16. Two different CMDs for our sample, both using the (J - [VTl]) colour, that separates well the different populations. Left: Absolute 
magnitude in [W4], dominated by the outermost layers of the stellar envelope. In this plot, symbol size is a function of LC (i.e larger symbols 
imply brighter classes). Right: Absolute magnitude in the J band, for which the bolometric correction is lower. 


tion of almost 3 magnitudes in the /-band, much larger than 
the typical variations for RSGs, A/ < 0.5 mag (Groenewe- 
Igen et al.|2009[ ), and in concordance with the severe spectral 
variations. 

In consequence, despite its high luminosity, M^s - 
-9.55 mag, more typical of RSGs than AGB stars, we have 
to conclude that HV 838 is a very luminous AGB star. This 
is in agreement with its position in the period/luminosity di¬ 
agram ([ Wood et al. II 1983 1) and the presence of a significant 
Li I 6707A line ( [Smith et al.|l9^ . 

HV 1956 ([M2002] 5873 8): This star i s a well know n long- 
period Cepheid variable ( |Butler|[ 
very long period of - 


;r||T97^ 

210 days ( [Soszynski et al.|20 


Eggen|1977 


Thi s 

star has been classified as MO-11 by (jPrevot et al.||1983|l, 
G2 Ib by |Wallerstein| ( |T984) l and K21 by |Massey|P0021 l. We 


with a 


have observed it in all three epochs, finding spectral types 
G5Ia in 2010, G7Ia in 2011 and G4Ia in 2012. This is a 
very peculiar object and while we cannot discard the possi¬ 
bility of very large spectral variations, we can conclude that 
it spends most of its time in the G spectral type, as expected 
for a Cepheid variable. 


In addition, our red spectra show inverse P-Cygni emission 
for most of the lines in 2010, no emission in 2011, and 
Cepheid-like profiles for the Ca triplet lines only in 2012. 
Finally, we have to note this star has a(J - H) ~ 0.4. There¬ 
fore, in Fig. it is the only G star with the expected values 
for giants (the solid red line). 

- SMC091 (PMMR 10): This star shows a low radial velocity 
(fhei = 54.5 kms^'), a value well under our velocity thresh¬ 
old to be considered a member of the SMC. However, our 
classification for it is K4Ia. Its high luminosity cannot be 
doubted and therefore it is a RSG, as first suggested by |Pre- 
|vot et al.| ( |1983 )l. It has Ks - 9.062 mag, typical of SMC 
RSGs (see Fig.[^. Therefore, despite its anomalous velocity, 
it is an RSG at the distance of the SMC. We speculate that 
is may be a runaway, ejected from a cluster or, more likely, 
from a binary in a supernova explosion. 

- SMC121: This star presents very atypical MIR photometry, 
with low photometric errors. It has a M[wa] - -15.7 and 
([1T3] - [1T4])= 4.4, indicative of an extreme dust enve¬ 
lope. Both values are close to being the highest among all 
our RSGs and AGB stars. The spectrum of this star presents 
no anomalies, and we have classified it as G8 Ib. Inspection 
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of the WISE W4 images for the area show that the star is pro¬ 
jected on a patch of extended bright emission, associated to 
the H II region LHA 115-N 23. The anomalous MIR colours 
are undoubtedly due to contamination by the extended emis¬ 
sion. The nearby RSG SMC 123 is also likely contaminated 
by emission from this Hii region and/or LHA 115-N22. For 
this reason, both stars have been removed from the plots. 
SMC145: Its NIR spectrum is similar to a Carbon Star, but 
the Ca triplet is still visible. Its optical spectrum is not that 
of a Carbon star either, but it does not look similar to any of 
our standard or reference stars. We speculate that it may be 
related to S stars. 

SMC169: This very late star displays an optical spectrum 
around M8.5 III. However, it has a Mkj = -10.2 mag, more 
typical of a RSG than a giant star. Soszynski et al. ( |2009| l 
give a main period of 1062 days, with an amplitude in the 
I band of 2.7 mag, much larger than the typical values for 
RSGs, (~ 0.5 mag; Groenewegen et al.|2009| l. We have also 
examined the NIR spectrum, which shows many metallic 
lines despite its late spectral type. Therefore, we consider 
it as very luminous AGB star, and we have classified it as 
M8II. 

SMC283 (CM Tuc); This is the brightest star in our sample, 
with Ks = 4.785 mag. It also has the only negative radial 
velocity and the latest SpT (M6 lab) in our sample. Attending 
to its velocity and brightness, it cannot belong to the SMC. In 
fact, it was previously identified as a foreground star because 
of its velocity by |Prevot et al.| ( |T983| l. 

In mid or late-M luminous stars, the rise of TiO bands erodes 
the continuum, weakening or erasing the atomic lines. In 
consequence, the exact luminosity class of this object can¬ 
not be ascertained. Being a foreground star, we have no in¬ 
formation about its intrinsic luminosity. Therefore, even if 
morphologically it shows the properties of a supergiant, it 
may be an AGB star. 

The average intrinsic Ks magnitude for other stars with 
similar spectral type (M51) in our SMC sample is Mks = 
-10.2 mag. If we consider this as the intrinsic magnitude for 
our star, we obtain for it a distance modulus of 14.5 mag, 
i.e. 10 kpc, too far from the SMC to have any relation with 
it. However, as a galactic star, the reason for this location and 
velocity remains without explanation. 

SMC311 (HV 12149): This very late star shows a spectrum 
about M8III. However, it has a Mkj = -10.4 mag, more 
typical of RSGs than giant stars. Soszynski et al. ( |2009| l give 
a main period of 769 days, with an amplitude of 2.3 mag, 
much larger than the typical variations for RSGs. We have 
also examined the NIR spectrum, which shows many metal¬ 
lic lines despite its late spectral type. Therefore, we consider 
it as very luminous AGB star, and we have classified it as 
M8.5 II. 


- SMC401 (HV 2112): With a special type about M5.5 II, this 
object has Mk^ = -10.3 mag, again very bright for a giant. 
However, OGLE has classihed it as an unresolved multiple 
star, perhaps explaining its atypically high brightness in the 
Ks band. In any case, this has to be a very luminous AGB 
star. 

- LMC039: Its velocity is higher than the assumed threshold 
for the LMC. Its spectrum, however, corresponds without 
a doubt to a RSG, while its brightness is inside the typical 
range for RSGs in the LMC. In consequence, we consider it 
as an LMC RSG with peculiar velocity. 

- LMC074 (HV 2572): This object was proposed by |Wood~& 
|Bessell|([l985|l as a candidate low-luminosity RSG, because 
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of a photometric period of only 201 days. The spectroscopic 
data disprove this possibility. Our spectral type, combining 
the blue and infrared spectra, is M7 III, while [Smith et~ar 
(1995 1 report the detection of a strong Li i 6707A line. Even 
though Groenewegen et ah] ( 2009} report a photometric pe¬ 
riod of 312 days, OGLE III data ( [Soszynski et al.||2009| l 
give a main period of 605 days, more in line with the lumi¬ 


nosity/period relation for AGB stars. Also, Soszynski et al. 


(2009 1 observed an amplitude in I of 2.5 mag, too large for a 
RSG ( Groenewegen et al.|2009| l. HV 2572 is therefore a lu¬ 
minous AGB star, with Mkj = -10.1 mag (above the upper 
Mks limit in the period/luminosity diagram of Wood et al. 


[1983[). Therefore we have opted to classify it as M7.511-111. 
LMC169 (HV 2670): This star has a radial velocity below 
the threshold adopted, but its spectrum shows clear RSG fea¬ 
tures. Its brightness is also typical for a LMC RSG. There¬ 
fore, we consider that this star is a LMC RSG with peculiar 
velocity. 


We have also found emission in Baimer H lines for some 
of our RSGs. As they present no other typical nebular emission 
lines, we consider this emission as intrinsic. On the other side, in 
the IR spectra there are no emission lines or any other peculiari¬ 
ties. In the SMC, stars with Balmer line emission are: SMC374, 
SMC372, [S84d] 105-7, and [M2002] SMC 8324, 8930, 9766, 
13472, 18592, 23463 & 55355. In the LMC sample, we find: 
LMC122 and [M2002] LMC 143035 & 148381. 

A number of targets show evidence for a blue companion in 
their optical spectra. We list: 

- [M2002] 55933: The blue end of its optical spectrum is dom¬ 
inated by the signal of an early B star. We have checked avail¬ 
able images in the U band and found a large number of blue 
objects around our star. Therefore, this early-B star might be 
a visual companion. 

- [M2002] 67554: For wavelenghts shorter than 4300A the 
spectrum is dominated by the flux of an early-B star. In the 
available image in the U band, there is a blue star 4" away. 
The possibility that some of the flux is collected by the fibre 
cannot be discarded, and so the blue star might be the visual 
companion. 

- [M2002] 51906: This star has an H5 absorption line stronger 
than usual for its SpT. It may be caused by a physical early- 
type companion, as no blue stars close to this RSG are seen 
in the l/-band image. 

- YSGOlO: The blue end of its optical spectrum is dominated 
by the flux of a B star. Given the shape of its lines, it has to 
be a fast rotator. As in the t/-band image there are not other 
blue stars close to it, the B star may be a binary companion. 
This object presents emission in the Balmer lines, but there 
are no nebular emission lines, nor emission features in the 
nIR spectrum. 

- [M2002] 169142, 168047 and LMC238: In all cases, the 
blue end of the optical spectrum is dominated by the flux 
of an early B star. As these stars are in the middle of clusters 
(H88 298, KMK88 91 and BSDL 2654, respectively) the B 
star is probably another spatially-separated member. 

- LMC172: The blue end of its optic spectrum shows the in¬ 
print of an early star. The l/-band image does not provide 
definite information, and this star may be a visual compan¬ 
ion. As the S/N is low for this possible companion, a more 
detailed classification is not possible. 

- LMC049: This star is in the middle of the cluster NGC 1967, 
and the blue end of its optical spectrum is dominated by the 
flux of a BO-B1 star. 
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- LMC239: This star is in the middle of the cluster H88 301, 
and its blue spectrum shows traces of an early B star. 

- LMC256; This star is in the middle of the cluster H88 308, 
and its blue optic spectrum is dominated by a ~B 1 III star. 

- SMC099; An early-B star appears in its blue spectrum. How¬ 
ever, in the U band image there are no bright sources close to 
this star. Therefore, this early-B contaminant may be a phys¬ 
ical companion. 

- LMC062: This star is in the cluster NGC 1983. Its optical 
spectrum is completely dominated by an B91 star, but its red 
spectrum is a blend of a bright K-type supergiant and an ear¬ 
lier component, probably a bright G o F-type star. 

- LMCllO: Despite the late spectral type of this star (M5Ib- 
II), it presents strong Balmer lines in the blue sectrum. Since 
no He I lines are seen, this early-type companion may be a 
late-B or early-A star. In this case, its LC should be at least 
II to be observable in the LMC. There are no indications of 
a cluster close to this star. 


5. Conclusions 

We have performed a pilot study in Large and Small Magellanic 
Clouds, aimed at their red supergiant population. Over a set of 
photometrically selected candidates, we have performed a de¬ 
tailed spectroscopic analysis, deriving spectral types, luminosity 
classes and line-of-sight velocities for all the observed targets. 
Once classified and with the available photometry, we show that: 

- There is a large population of supergiants in both clouds, 
largely in the dim end of their brightness range, that remains 
to be observed. 

- There is no purely photometric criterion capable of sepa¬ 
rating completely different populations, and when selecting 
RSGs, we will always have to choose between the complete¬ 
ness of the photometric sample and its cleanliness. Due to 
the fact that instead of clear-cut borders, the transitions be¬ 
tween different populations are gradual, there will always be 
some AGB and carbon stars that will appear as interlopers. 

- It is possible, nonetheless, to use the synergies between near 
and mid infrared photometry to open avenues to much more 
efficient selection criteria. 

- The completeness of the photometric selection criteria is a 
function of spectral type, and in particular there is a loss 
of efficiency for red supergiants with the earliest and latest 
types. This has to be weighed in whenever drawing conclu¬ 
sions about their relatives abundances in several astrophysi- 
cal contexts. 

- Mass loss becomes important only for supergiants later than 
K5, although it is not ubiquitous and at each spectral type 
there will stars with no or very little apparent mass loss. 

- The thermal behaviour of the dust inhabiting these expelled 
outer layers seems to be similar for all the supergiants, inde¬ 
pendent of spectral type and luminosity class. 
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